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Abstract 

A  aa tbs  as t leal  expression  is  written  Seen  for  the  changes  in  transient 
frequency  and  daeplng  which  result  whan  the  feedback  branch  of  a  single 
loop  syatea  is  altered  slightly.  The  result  is  applied  in  a  general  dis¬ 
cussion  of  the  accuracy  of  aerodynamic  simulators  with  numerical  examples  of 
roll  el aula t lon^  and  suggested  experiments  to  evaluate  the  approximations  in¬ 
volved  in  practice, 

^  Stesaaafi. 

A  simulator,  when  used  to  test  a  missile  guiding  control,  completes 
e  feedback  loop  by  supplying  the  ^  (feedback)  part  of  a  lAi-Beta  (y44  -  @  ) 
system,  where  represents  the  amplification  in  the  forward  transmission  path. 
The  transient  response  of  such  a  system  to  any  disturbing  influence  is  a  sat  of 
oscillation#  2a,  at  coaplax  frequencies  +  j  ,  wbart 

are  the  root*  of  tha  expression  (  1  -M  §  ),  i-,  e.  solutions  to  tbe  equation 
M  p  m  •  +  j  o  0  then  all  valuas  of  X i  are  negative,  the  syatea  la  said 

to  be  stable,  in  that  transiant  oscillations  will  decay  with  time  after  they  are 


// 


-2 


initiated.  The  problem  ir.  control  design  is  to  make  all  the  sufficiently 

negative  so  that  the  transient  component  os  dilations  will  decay  rapidly,  and 
quiescence  will  be  restored  sufficiently  soon  after  any  disturbance. 

Special  consideration  rust  be  given  to  the  effeot  of  a  simulator  which 
does  not  provide  exactly  the  same  as  the  dynamical  system  it  simulates,  but 
only  approximately  so.  The  purpose  of  simulation  experiments  is  to  determine 
the  important  -'s  and  CJ3  '*  of  the,/**  ^3  system  which  will  obtain  when  the 

missile  is  in  flight,  and  it  will  be  shewn  below  that  the  precision  obtained  for 
and  60^  on  a  percentage  basis  is  not  at  all  the  same  as  that  obtained 
for  the  absolute  value,  j j  ,  on  a  percentage  basis.  An  expression  will  be 
developed  for  snail  deviations  from  desired  flight  values  in  terms  of  ft 
variations  from  actual  flight  behavior. 

Let  9e  ■*  if  +  J  *|  be  s  small  perturbation  in  ft  so  that 
equals  feedback  given  by  simulator,  ft0  equal*  feed  bade  required  by  the  equation 
of  moti on .  Then  £•<».-  h  and  observed  transient  frequenoiee  will  be 
found  fer  ft  $  ~  &  or  :>r  f  ~/*  ft o*  S’ ft l  iMt«d  of  sere. 

The  values  corresponding  to  ^  t  —yU  ^ 


The  values 


can  be  found  o-  observing  that 
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as  long  as  is  not  a  singular  point  of  the  !"/** 
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•ufficiently  small 


•  fih 


sufficiently  small.  Strictly 


•pocking,  ff  ahould  bo  evaluated  at  tho  observed  "frequency"  #■  in  thi« 

P*  '  /?  J 

aquation,  but  it  is  doubtful  if  tho  refinement  of  distinguishing  uetwoau  f  ej 
,Dd  R  '  will  bo  worth  whilo  in  view  of  tho  othor  unoor taintioo  in  tho 

r* 

•ituation,  thoroforo  thi«  distinction  ia  ignored  henceforth. 

An  example  of  tho  foregoing  can  bo  found  in  tho  desi&i  and  uso  of  a  roll 
•  inula tor.  the  equation  of  potior,  of  a  miasilo  in  roll  ia 

A<p  +  G-  <p  *  -  //<< 


whore 


oC 


•  angle  of  roll 

•  angle  of  control  fin 

•  rotational  inertia 
"  damping  coefficient 

•  fin  effectiveness  coefficient 

The  (2  ia  ■~jr'  computed  from  thia  equation,  and  the  y**  ia  — ^ 
computed  from  the  internal  oortrol  system.  Tha  transient  osoillationa  are  given  by 
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Gr 

H 
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and  /<f  ~  ( -fj 


*f+*r 

thu.  ±J  efty  for  an  idealized  control 

where  f*'  ■  conatant  iiepondent  of  ^  . 


Two  examples  will  bo  treated  quantitatively  for  roll  simulation  i 


06.no*  for  roll  aimlation 


(tbo  idealised  oom,  uaaful  for  eiaolator  oalibratlon) 
nlotloa  p  then 


r  '  (A*1"-*-  S"f») 


+  &) 


-  Gr  .  -  vl  /*H  ^  /  A-  \** 
evaluated  at  y/a'^  =  j,  —  J  f  '"a”'  i  >  A  / 


Tfao  oboerved  damping  oonatant,  y\,  f  la  accordingly  raduood  la  tto  ratio  '  oJ~  Q 
and  tbo  observed  frequency y  Ce>;  ,  of  trmnalaot  ooolllatlogo  bgr  |  - 

*  A 

free  tbelr  accurate  values ,  duo  to  tbs  approximation  ^  -f-  |  ^  ^  q 


For  tbo  following  nuaerical  values ; 
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S'  »  1/6 

^  •  yA 0  lboino*  *  8.93  dynaalc  unit* 
H  -  36900  lb* in. /radian 
£  •  4)1  lb^iiio/rad  /aaoo 


it  la  found  that 


=  f^r-ihr  -  ,3+rs~  ' 

4*  ”  g*  -  -22.?  sec.  * 
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=  Co*  (/-  .3|e 

Showing  that  far  X^  and  to/ to  ba  oorract  within  5#  wa  aust  bmrm/f/ <,510  =  .  iLi 


/f/< 


Si  no* 


^.5«  tj  u  j,  tha  peroantago 
'4^T°bat  thia  appliaa  for  y**  *  1/6  only, 


oner  in  /£/  want  ba  laaa  than  0r4»- ■  ^ftr°bnt  thia  appliaa  f oar/44  •  1/6  only 

6 

far  daoping  laaa  than  about  *706  orltieal,,  tha  parcantaga  arror  raaulting  fraa  tha 
aaa  of  an  approxlaat*  alaulatar  variao  with  9  ao  that  for ■  1/2  (a  rain* 
vhlah  baa  baan  auggaatad)  tha  praaialon  la  j@i  auat  ba  thrao  tiaaa  aa  good  aa 

ab ora,  in  apita  of  tha  aaallar  valna  of  ft]  whiah  obtain*  for  tha  largar  0 

r<- 

A  portion  of  tha/*  |?  plana  for  thla  axaapla  la  plotted  in  Figaro  1,  with 
tha  line*  at  CO  •  13.43  andX*  -42.3  ■  pool  ally  dark#  nod  Tha  branab  paint  at/*4  /  * 

1.36  ♦jo  la  ooanaatod  with  tha  oxitioal  daaplng  eharaetariatiaa  of  tha  ay  a  tea, 
and  oaanot  ba  loeatad  pradaoly  without  knowing  tha  analytic  fora  ofy**  ,4  aa  a 
fanatloa  of  .,ja  • 


University  of  Virginia  of  Internal  Control 
Reference  (a)  gives  an  analytic  expression  for  tne  dynamic  behavior 
of  the  University  of  Virginia  servo  as  a  third  degree  equation  in  .  This 

squat:  n  is  obtained  oy  analyzing  tne  dynamics  of  the  control  motor,  but  ignoring 
inertia  l«g  in  the  hydraulic  system.  While  the  later  data  of  Reference  (b)  can 
be  fitted  much  better  by  a  foirth  degiev  equation  in  ^  (vnich  includes  a  saall 
correction  term  for  hydraulic  fluid  inertia i,  the  cubic  expression  of  Reference  (a) 
will  be  used  in  this  illustration- 

Thus,  for  a  control  whose  internal  behavior  is  specified  by 

y*:  y*- _ __ 

-fJ+ ,v v -pi‘+  -y-  %». o i.  x  io * 

(See  Reference  (a)  pp.  and  2*t) 


and  whose  external  feedback  is  given  by 


B  - - ftu 
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(Sane  aa  Oase  I) 


we  have 


/i.  g  2-oT *  /c _ _ _ —  - - — 

if  tne  linkages  are  ao  proportioned  that,  for  very  alow  variations, 

This  numerical  scale  corresponds  to  that  of  Ca»e  I  above- 

Tne  transient  oscillation  "frequencies*  are  given  by  the  root*  of 


i  Ifltto 


lf*  +  ZO<\xio,*'fS+*Sc.:  '°S  f 


-i  i  li,  \  •  ’  *>  i 


These  roots  are 


-  s*  i  ■-  si  z  , 


and  ..  J.'  ;  I  - 


The  transient  term  of  moat  importance  la  the  oscillation  mode  at 

frequency  16, 1  =  2-5  epe,  since  not  only  are  the  higher  order  roota 
2 

cheracteriaed  by  suen  greater  damping,  but  the  amplitudes  of  the  oscillations  at 
frequencies  of  larger  are  much  less  than  that  of  the  lowaat  aode  in  gust- 


initiated  transients 


.7- 


Tbe  value  of 


<f  to 


UL 


■f 

evaluated  at  each  of  t»«e  five  root*  above 

-L-  !a 


Therefore  at  =-i.7+/L  <j  ,  ^  ~^-3Ho  4  7  7*1°  J 


Hence 


r 3  /: 


+i 


fo  -  <WJ 


775T,  »  ih  *  ~fc~  (■"+  ”j  ) 


'J 

-  V/3o 


at  ft  »  -L-7  41L lj ,  ~  jts-  -.  /aij 

•  '.  //=  -  |?€  +  °  St,  J  =  -4.7+  /<.  lj  -  $%  (M*  +  °  s<0p 


But  since 


^  Zf+J\‘  fc  ='<-7(,  +  ‘Z^^'  tfS3  J  (l-  *3*  7.-- 

and  the  percentage  error  involved  la 


in  the  duxt'.  ,.-;  l'actor  and 
lr.  frequ,  .**:> 


-J.i"  ^  7 

In  Figure  ^  la  plotted  a  portion  ci  tee/*  plane  ttt  region  around 

— 4. ^ f  j  It  is  not  consistent  witn  t.»e  title  of  this  paper  to  explore 

the  irregular;  ties  around  the  branch  points  &t,-3(,.S  4  J°,~  IO.4L+J0,  -lSMt±l7sj 
since  these  are  not  close  to  tne  "resonant"  points  given  by  ^  a  /+j  « 


B<  Daalim  q-’  Hbcperljej;te_ 

The  aatheuetlcal  treataent  in  Section  1  of  this  report  has  been  carried 
somewhat  beyond  the  margin  of  utility  for  the  case  of  the  University  of  Virginia 
control,  In  vie*  of  the  iaoerfect  approximation  to  Its  dynamic  behavior  fur¬ 
nished  by  a  characteristic  equation  of  degree  no  higher  then  the  third,  It  la 
the  purpose  of  this  section  to  suggest  experiments  which  will  give  e 


reasonable  basis  for  eetiaating  tbs  ^  and  /j»c  appropriate  to  a  given 
experlaental  setup  for  any  control  of  the  general  types  heretofore  suggested., 

In  flection  A,  the  observed  value  of  transient  ooaplex  frequency  , 

the  accurate  value  j°L  ,  and  the  slwolator  error  f  t  are  shewn  to  be 
related  by 

,  /‘lb 
A  *  a  *  irnf^lT 

‘  °  lon«  ■*  and  !@*/ 
are  ssall  enough  so  that  the  first  two  tens  of  the  Taylor  series  fora  a 
sufficient  approximation  to  the y'**'  p  function  In  the  neighborhood  of  the 
Hyqulst  point.  To  the  sane  degree  of  a pproxi nation  the  difference  oan  be 
written 


1  -zrrxf 


li 


Is  thus  given  In  terns  of  experlaentally 


M)L 

•  r* 

end  It  Is  to  be  observed  thet 
deternlm..  quantities, 

i 

The  value  of  yj  sen  be  observed  directly  by  measuring  frequancy 
and  logarl thalc  dacreaant  of  transient  oscillations  initiated  by  e  step- 
function  signal  fad  into  the  system.  By  vatylng  the  nuasrlcal  coupling  factor 
between  the  control  and  the  (Isolator,  thara  rasulte  a  change  in  aeale  of  the 
/?,  plot,  so  that  the  affect  Is  one  of  at  -mg  the  By quiet  point.  If  the 
numerical  aagnltuda  of  ^  is  thus  changed  by  the  ratio  ~ft  ,  the  resulting 
corresponds  to  that  obtained  from  the  original /**  &  plot  at  the  point 
j  °  hence  the  observed  ^  'fi  corresponds  to  a  J.  of 


i 


'  -  It 


The  » a. 


ta  seat  Dc  determined  by  coepuuing  (mm  the 


0  ' 

aotior.  equations.;  U.e  |  appropriate  to  the  observed  ,  arid  coapailng  It  to 

the  observed  aioslle  aoveaent  ratio  from  the  oscillegraas  of  ti*e  siMuiator- 

control  surface  aoticn 


n'.trol  system  The  value  of 


1  _l — 7 

'  p  will  be  given  by  J 


It  is  appropriate  to  reeark  that  the  values  of  frequency  and  logariti»lc 
decrement  ere  not.  very  easy  to  measure  accurately  on  the  oscillograms  in  tha 
neighborhood  of  critical  daaping  It  will  usually  happen,  however,  that  consider 

i 

able  errors  in  ^  —  -p  can  be  tolerated  in  rne  neighbornoou  of  or: heal  damping 
just  because  tha  cnange  thereby  introduced  in  the  aeawuralla  benevior  o?  the  systaa 
will  not  be  important^ 

A  eore  fundamental  objection  is  that  the  aatneeatical  treateer.t  of  this 
paper  is  based  on  the  assumption  thaty^  exists  as  defined,  and  is  an  ar*i;/\ie 
function  of  Tne  usual  reaarke  aoout  linear  ap:-ro:  ‘  V.-ne  to  non-linear 

systaas  can  be  aade,  however,  since  tiie  anc.jysia  refers  only  to  a  differential 
region  around  tha  NyquJst  point,  hence  is  an  appropriate  epproxi nation  for  any 
ta*  which  is  continuous  enough  in  its  dynamic  behavior  to  warrant  analysis 
fcy  differential  means  at  all.  The  extant  wr.ioh  ?  depa  rts  from  being  analytic 
can  probably  b*  explored  by  introducing  phase  shift  as  well  as  magnitude  shift 
in  cnang  ng  the  scala  of  theyA  8  plot  (see  page  ),  tnus  evaluating  ~ ) 

:  o-  other  than  real  values  of  d*  • 
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